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I.  IHTBOD OCTION 


accurate  and  reliable  numerical  weather  prediction  is  a 
major  goal  fee  civilian  and  military  atmospheric  studies. 
Verification  studies  are  an  important  aspect  of  numerical 
model  develcpaent  as  they  isolate  model  strengths  and  weax- 
nesses  and  identify  errors  which  must  be  eliminated  tc 
extend  the  present  limits  cf  forecast  shill. 

This  thesis  will  investigate  the  Naval  Operational 
Global  Atmospheric  Prediction  System  (NOGAPS)  model's  abil¬ 
ity  tc  forecast  the  pcsiticn  and  movement  of  500mb  waves  ir. 
the  mid  la tit  ides  of  the  Northern  Hemisphere.  NOS  APS ,  which 
became  operational  in  September,  1982  at  the  Fleet  Numerical 
Oceanography  tenter  (FNOC)  has  been  a  project  of  the  Naval 
Environmental  Prediction  Eesearch  Facility  (NEPHF)  since 
1976.  The  long  term  goal  of  NOGAPS  is  to  achieve  a  medium 
range  (7-  (C  days)  numerical  forecast  capability  for  the  U.S. 
Navy  (Bcsacnd ,  1981)  .  Verification  cf  NOGAPS  is  very  impor¬ 
tant  at  this  stace  due  tc  the  research  type  background  cf 
the  mcdel .  The  NOGAFS  forecast  model  is  a  modified  form  of 
the  CCLA  general  circulation  model.  Consequently,  it  has 
not  had  the  benefit  of  aany  prior  intensive  verification 
studies  as  have  most  second  and  third  generation  operational 


Becent  studies  have  raised  serious  questions  concerning 
the  relative  accuracy  of  the  planetary-scale  (zonal  wavenum¬ 
bers  1-2)  wave  action  forecasts  at  SGOmb  produced  by  numeri¬ 
cal  weather  prediction  models.  Lambert  and  fierilees  (1978) 
concluded  the  planetary-scale  motion  predictions  of  a  spec¬ 
tral  numerical  model  were  poorer  with  respect  :c  persistence 
than  long  (wavenumbers  3-5)  or  medium-scale  (wavenumbers 
6-10)  moticns.  Miyakcda  et  al.  (1972)  examined  the  GFDL 
model  performance  over  a  series  of  winter  forecasts  and 
found  similar  error  characteristics  in  the  planetary-scale 
forecasts.  This:  study  will  verify  NOG  APS  planetary  wave 
forecasts  ar.c  cti.er  wave  regimes,  large  ar.d  medium-scale. 

50Cmb  data  from  both  the  forecast  and  the  analysis  are 
spectrally  deccmpcsad  intc  its  basic  components  of  amplitude 
and  phase  then  organized  into  planetary,  long  and  medium- 
scale  waves  for  separate  analysis  cf  therr  atmospheric 
behavioral  characteristics.  If  one  assumes  spatial  coupling 
(the  interacting  of  different  scales)  ,  then  model  weaknesses 
in  forecasting  'long'  wave  motion  will  be  reflected  in  the 
accuracy  of  'short'  wave  notion  forecasts.  The  purpose  of 
this  study  is  tc  perform  a  thorough  analysis  cf  the  fore¬ 
casts  of  each  wave  grouping.  Somerville  (1980)  examined  the 


planetary-scale  wave  forecasts  cf  a  primitive  equation 
numerical  model  in  bcth  a  global  and  hemispheric  configura¬ 
tion.  He  concluded  that  the  global  version  was  markedly 
mors  skillful  than  the  hemispheric  version  especially  in  the 
latter  part  cf  the  five-day  forecast  period.  While  the 
present  study  will  consider  only  a  global  version  cf  NOGAES 
the  results  will  lend  themselves  to  an  application  of  Somer¬ 
ville's  work  tc  NOGAPS. 

Errors  in  phase  speeds,  amplitudes,  development  of  ba re¬ 
el  ir. ic  systems,  damping  of  smaller  scale  features  and  plane¬ 
tary  wave  structure  become  clearer  when  the  forecasted 
variables  are  analyzed  over  time  via  a  Hovmoller  or  trough- 
ridge  plot.  5=  umhfefr. er  and  Downey  (1978)  u-ilized  this 
technique  tc  inter compare  the  forecasting  skills  of  three 
numerical  weather  prediction  models.  They  constructed  lon¬ 
gitude-time  plots  along  a  fixed  latitude  to  better  under¬ 
stand  the  time  evolution  cf  the  forecast  fields.  Forecasts 
which  originally  appeared  quite  good  were  revealed  tc  have 


amplification  errors 

in 

many  of  the  transient 

systems.  A 

spectral  breakdown 

cf 

the  Hovmoller  plot 

in  planetary 

(wavenumbers  1-3), 

large  (wavenumbers  5-7) 

and  medium 

scale-waves  (wavenumbers 

8-12)  were  utilized  ir 

i  their  study. 
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One  has  only  to  consider  the  complexity  of  a  typical  500mb 
wave  pattern  tc  appreciate  the  benefits  of  this  approach. 

The  specific  objectives  of  this  thesis  are: 

(1)  Development  of  a  spectral  verification  program  using 
wavenumber  grouping  and  Hcvmoller  plots  for  tropospheric 
height  data, 

(2)  Application  of  the  prograa  to  interesting  NOGAPS  fore¬ 
cast  situations  at  530mb, 

(3)  Use  of  wavenuafcer  analysis  to  verify  NOGAPS  forecasts 
of  planetary,  long  and  medium-scale  waves. 

Chapter  2  reviews  the  main  features  of  the  NOGAPS  model. 
The  pattern  decomposition  program  is  described  and  demon¬ 
strated  in  chapter  1.  Chapter  4  presents  the  basic  Charac¬ 
teristics  of  the  case  study  periods  and  their  subsequent 
atmospheric  behavior.  finally,  chapter  5  summarizes  the 
results  of  the  forecast  ccnpariscn  and  makes  suggestions  for 
further  research. 
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II.  NCGAPS  SUMMARY 

The  Naval  Cperaticn<il  Global  Atmospheric  Prediction  Sys¬ 
tem  (NOGAFS)  was  developed  by  the  Naval  Environmental  Pre¬ 
diction  Research  Facility  (NEPRF)  to  be  run  at  the  Fleet 
Numerical  Oceanography  Oerter  (FNCC).  The  objectives  cf  the 
NOGAFS  project  are  to  improve  tropical  and  Southern  Hemis- 
peric  forecasting  and  achieve  aid-range  (5  -  13  days)  global 
numerical  forecast  capability  to  support  3.S.  Naval  opera¬ 
tions  (Rosmcnd,1981)  . 

The  objective  analysis  method  used  by  NOGAPS  is  a  varia¬ 
tion  cf  the  successive  corrections  technique  (Barnes,  1964) 
which  discriminates  among  observing  systems  according  to 
their  accuracy.  While  the  final  analysis  is  not  as  sophis¬ 
ticated  as  the  optimal  interpolation  method  it  requires  far 
less  computer  resources. 

Glooal  mass  and  wind  fields  are  analyzed  independently 
with  optional  data  bogusing.  Eata  rejection  criteria  are 
intentionally  lcose  since  the  initialization  scheme 
(Barker,  1 981)  removes  the  gravity  waves  excited  by  bad  data. 
This  technigue  prevents  inadvertent  rejection  of  the  good 
data  associated  with  rapidly  developing  systems. 
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A  calculus  cf  variation  method  combines  the  results  of 
the  wind  and  aass  analyses.  The  hydrostatic  equation  per- 
lits  substitution  of  teiperature  for  geopotential  which 
gives  the  forecast  model  three  diaensicnality .  The  fora  of 
HOGAES  initialization  used  in  our  study  produces  balanced 
geopctencials  (temperatures)  and  non-di  verg ent  winds  which 
are  interpolated  directly  to  the  medal’s  sigma  coordinate 
system.  The  method  which  became  operational  in  December, 
1982,  defines  gaopotential  and  wind  as  objective  analyses 
correction  fields.  The  resulting  balanced  corrections  are 
interpolated  tc  the  sigma  surfaces,  as  in  the  first  method, 
and  added  tc  the  first  guess  forecast.  This  method  pre¬ 
serves  the  model-generated  divergence  and  minimizes  the  ver¬ 
tical  int erpclatica  error. 

The  first  guess  forecast  of  specific  humidity  and  the 
prognostic  planetary  boundary  layer  (?BL)  variables  are  used 
as  initial  values  for  the  subsequent  forecast  cycle.  The 
six- hour  update  cycle  is  preferred  over  the  twelve-hour  con¬ 
figuration  since  less  tine  interpolation  of  asynoptic  data 
is  required.  However,  the  shorter  cycle  is  more  subject  to 
initialization  shcck  and  relies  heavily  on  the  initializa¬ 
tion  technique  to  control  the  noise.  Failure  to  do  sc  may 
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result  in  subsequent  first  guess  errors  sufficiently  large 
to  cause  rejection  of  good  observational  data. 

lhe  NOGAES  model  is  a  modified  version  of  the  UCLA  gen¬ 
eral  circulation  model  (GCM)  described  by  Arakawa  and  Lamb 
(1977).  lhe  2.4  by  3.0  degree  horizontal  grid  uses  scheme  C 
staggering  cf  the  variables  and  conserves  energy  and  enstro- 
phy  when  the  flow  is  roniivergent.  All  dependent  variables 
except  vertical  velocity  are  defined  in  the  middle  of  each 
of  the  six  vertical  sigma  levels  (Fig.  1)  . 

Moisture,  mcist  static  energy  and  momentum  are 
assumed  tc  be  well-mixed  in  the  PSL.  The  equations  govern¬ 
ing  the  diahatic  processes  ate  based  cn  the  method  cf  Dear- 
dcrff  (1972).  These  equations  represent  the  conservation  cf 
mass  in  the  PEL,  temperature  and  specific  humidity  inversion 
strength  and  the  jump  in  momentum  across  the  inversion.  A 
special  feature  of  the  mcdel  is  the  variaole  depth  of  the 
PBL  ir.  the  lowest  layer  which  has  a  capping  inversion  cf  a 
porocs  material  surface.  This  allows  entrained  mass  *o 
deepen  the  PBL  and  removed  mass  (by  cumulus  mass  flux)  to 
decrease  the  EEL,  The  fluxes  at  tha  top  of  the  PBL  ars  net 
a  constant  fraction  of  the  surface  fluxes  but  rather  a  func¬ 
tion  of  the  entrainment.  Dissipating  stratus  causes  a 
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particularly  lazes  flux  across  the  inversion  for  stratus 
capped  PBL's. 

A  significant  modification  to  the  SCH  is  the  cor.  strain  - 
ment  cf  the  PEI  to  remain  ir.  the  bottom  sigma  layer.  This 
eliminates  serious  ccmputationa 1  problems  yet  has  net  cre¬ 
ated  any  noticeable  adverse  effects  on  model  forecasts. 

NCGAPS  uses  the  Araka wa-S chu fcert  (1974)  cumulus  parame¬ 
terization  described  by  lord  (1978)  .  Environmental  air 
below  the  cloud  top  is  entrained  while  air  at  the  cloud  top 
is  detrained.  Arakawa  noted  that  large  scale  forcing  gener¬ 
ates  conditional  stability  whereas  cumulus  convec-icn 
destroys  it.  This  concept,  called  ' quasi-equilibrium' ,  was 
combined  with  the  'cloud  work*  concept  (Arakawa  and  Schu¬ 
bert,  1974)  to  produce  a  cloud  base  mass  flux  unique  to  each 
cloud  type. 

The  GCt!  radiation  parameterization  requires  only  fifteen 
percent  of  the  model's  running  time.  The  radiation  transfer 
equation  (Katayama,  1 S74  and  Sc hlesinger ,  1976)  uses  bulk 
transmission  functions  for  discrete  pressure  layers.  Net 
radiative  flux  at  ground  level  is  a  function  of  incoming 
solar  and  longwave  radiation  and  surface  albedo. 
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NCGAPS  is  run  on  the  CYBER  205  which  has  high  speed  sca¬ 
lar  and  vector  pipeline  processors.  Presently,  rns  forecast 
model  has  a  wall  clock  execution  time  of  thirty  minutes  per 
forecast  day  cn  a  2.h  X  3.0  degree  grid.  Any  significant 
increase  in  the  model's  performance  would  require  redesign¬ 
ing  the  diatatic  processes  to  include  more  vector! zation  for 
the  entire  model. 
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III.  SPEC? BAI  PS COMPOSITION  ROUTINE 


The  complex  appearance  cf  a  typical  500  mb  height  wave 
pattern  often  conceals  the  sources  of  forecast  error.  The 
pattern's  complexity  is  due  to  the  superposition  cf  plane¬ 
tary,  long,  medium  and  synoptic  waves.  Every  pattern, 
regardlass  of  its  intricacy,  can  be  represented  as  the  sum 
of  a  series  of  pure  harmonic  waves.  &  spectral  decomposi¬ 
tion  routine  will  resclve  a  coirplex  wave  into  its 
components,  that  is,  the  amplitude  and  phase  angle  cf  each 
wave  number  1,  2,  3,  etc.  The  spectral  analysis  used  in 
this  study  was  developed  at  the  Naval  Postgraduate  School- by 
Prof.  R.  T.  iiilliams. 

The  50C  at  wave  pattern  along  42N  (rig.  2)  for  18  Febru¬ 
ary  will  fcs  spectrally  decomposed  tc  illustrate  the 
technique. 

Tc  facilitate  this  illustration  and  the  foregcing 
sequence  of  illustrations,  Fig.  3  has  an  abcissa/ordinate 
scale  ratic  cf  25000km/6QCm  or  approximately  42000  to  1. 

Note  the  negative  longitudes  are  west  and  the  positive 
longitudes  are  east.  A  careful  examination  cf  the  heights 
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figure  2.  The  Height  contours  of  the  500mb  Surface  at  OOz 
18  Feb.  Contour  Intervals  are  120  Meters. 
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Figure  3. 


along  42N  suggest  a  relatively  significant  amplitude  of 
wavenumber  3  as  evidenced  by  the  prominent  ridge  locations 
in  the  vicinity  cf  120S,  0  and  90E.  Also  the  preponderance 
of  ridges  and  troughs  in  the  4000km  (50-60  degrees)  range 
are  the  signatures  of  wavenumbers  6  and  7 .  Since  a  minimum 
of  sii  points  is  required  to  adequately  resolve  a  wave,  the 
maximum  wavenumber  accurately  resclveable  with  a  72-point 
series  (five  degrees  of  longitude)  is  12.  A  horizontal  res¬ 
olution  of  three  degrees  cf  longitude  can  accurately  depict 
wavenumbers  up  tc  20.  Amplitude  estimates  are  difficult  ar.d 
finer  resolution  would  be  speculative  at  best. 

Ifce  numerical  results  of  the  decomposition  routine  are 
presented  in  Table  1  (  although  only  wavenumbers  1-20  can 
be  accurately  resolved,  the  tables  and  figures  in  this  study 
will  represent  wavenumbers  1  -  28)  .  The  phase  angle  ar.d 
amplitude  (in  meters  about  the  mean  height)  of  the 
wavenumber  1  through  28  components  are  shown. 

first  the  mean  height  was  calculated,  then  the  amplitude 
and  phase  angle  of  each  wavenumber  1  -  28.  The  amplitudes 
are  given  as  difference  values  about  the  mean  height  and  the 
phase  angles  are  in  reference  to  the  center  of  the  longitu¬ 
dinal  domain  (0  degrees)  . 
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TA&LE  I 

The  numerical  results  cf  the  spectral  decomposition. 

Spectral  deccapositior  of  500mb  wave  at  42H.  Amplitudes  are  in 
meters  and  the  phases  are  in  degrees. 


ASIJHH|1|C  SEAS 


HAVE  NOH 

1 

in&MDE 

PHASE 

-88.94 

2 

66.67 

27.38 

3 

92.17 

-157.89 

4 

24.9  1 

88.19 

5 

26.40 

126.31 

6 

46.9  2 

-41.75 

7 

57.17 

129.36 

8 

13.86 

134.82 

9 

13.63 

-0.63 

10 

10.19 

-  115.47 

11 

21.45 

101.96 

12 

10.36 

135.11 

13 

4.26 

-41.04 

14 

1.20 

142.10 

15 

6.06 

111.40 

16 

4.96 

125.36 

17 

9.54 

77.27 

18 

4.70 

1.36 

IS 

4.54 

64.28 

20 

2.9  1 

-52.99 

21 

1.13 

-96.49 

22 

1.56 

27.35 

23 

0.46 

-14.74 

24 

1.30 

-146.30 

25 

0.42 

-114.46 

26 

1.17 

136.61 

27 

1.3  0 

41.12 

28 

0.6  5 

-28.40 

A  graph  of  wave  number  versus  amplitude  (Fig.  4)  further 
clarifies  the  high  relative  amplitudes  of  wavenumbers  3,  6 
and  7  proposed  by  our  initial  visual  analysis.  The  compara¬ 
tive  importance  cf  wavenumber  11  is  also  apparent. 
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Figure  4. 
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The  series  cf  graphs  in  Figs.  5a-e  elucidate  the  signif¬ 
icance  of  "the  planetary  component  of  t ha  500ab  wavs.  All  of 
ohese  figures  craph  the  heigh*,  versus  the  longitude  cf  the 
original  wave  (solid  line)  with  various  wavenuirber  combina¬ 
tions  (dashed  line)  superposed  for  comparison.  Figs.  5a,  b 
and  c  chart  wavenumbers  1,  2  and  3  respectively,  while  d  is 
the  summation  cf  wavenumbers  1,  2  and  3. 


5e)  Summation  of  wavenumber  1-29  components 


Note  hew  these  three  waves  begin  to  resemble  the  origi¬ 
nal  pattern.  It  is  important  to  realize  the  summation  of 
the  full  series  cf  harmonic  waves  would  reconstitute  the  500 
mb  wave  eacepting  any  high  frequency  variations  (greater 
than  wavenumber  28)  present  in  the  original  wave.  Fig.  5e 
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represents  the  sub  of  wavenumbers  1  through  28  and,  or.  this 
scale,  appears  tc  coincide  with  the  test  wave. 

Ibis  analysis  approach  has  proved  to  he  a  useful 
research  tccl  when  applied  to  forecast  verification  (Baum- 
hefner  and  Downey, 1  978;  Somerville,  1 980)  .  Spectral  decompo¬ 
sition  cf  a  forecast  and  its  verification  analysis  allow 
direct  comparison  of  the  amplitude  and  phase  speed  of  each 
wavenumber.  The  dimension  of  time  further  strengthens  this 
analysis  *echnigue  via  a  Hovmoller  diagram.  These  trough 
and  ridge  lcngitude-t iaa  plots  along  a  fixed  latitude  depict 
-he  time  evolution  of  the  forecast  variables.  Fig.  6  com¬ 
pares  the  anclitude  and  phase  speed  of  successive;  24,  48  and 
72  hcur  forecasts  for  wavenumber  3  with  the  observed  atmos¬ 
pheric  behavior. 

In  the  next  chapter  this  and  other  cases  will  be  ana¬ 
lyzed  in  this  manner.  The  harmonics  of  each  forecast  are 
grouped  intc  planetary  (1 -3) ,  long  (5-7)  and  medium  (8-12) 
waves  befcre  comparison  to  the  verification  analysis  via 
Hovmoller  diagrams.  Beadiiy  discernaole  errors  in  ampli¬ 
tude,  phase  speed,  development  of  barociinic  systems,  damp¬ 
ing  cf  smaller  scale  features  and  planetary  wave  structure 
are  just  a  few  of  the  advantages  of  verification  using  spec¬ 
tral  components. 
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IV.  SPECTRAL  VERIFICATION  CASE  STUDIES 

In  this  chapter  the  spectral  verification  presented  ir. 
chapter  3  Mill  he  applied  to  three  NOGAPS  forecasts,  one  for 
three  days  and  two  extending  for  five  days.  All  forecasts 
are  Bade  from  00Z  data  anc  are  verified  at  24,  43,  72,  96 

and  120  hour  narks  of  the  forecast. 

4.1  Case  1  -  February  18  to  21,  1982: 

This  period  was  selected  cn  the  basis  of  its  rror.c  utcad 
planetary  and  long  wave  activi-y.  This  prediction  was  nade 
with  ar.  earlier  'coarse  resolution'  (4  degrees  let.  by  5 
degrees  long.)  version  of  the  NOGAPS  forecast  model.  The 
spectral  decomposition  was  applied  tc  features  along  42N 
averaged  between  38N  and  46N. 

The  initial  500mfc  pattern  for  case  1  was  discussed  in 
Chapter  3  and  is  shown  in  Fig.  2.  Prominent  features  along 
42N  include  planetary  wave  troughs  at  160»  and  45w  ar.d 
strong  ridges  at  170E,  HOW  and  15S,  suggesting  a  strong 
wavenumber  3  component  auplitude.  The  preponderance  cf 
ridges  and  troughs  with  wavelengths  on  the  order  of  50  tc  60 


degrees  are  the  signatures  of  wavenumbers  6  and  7.  E  s  c  a  1 1 
■chat  spectral  decomposition  of  the  heights  along  42N  con¬ 
fined  the  relative  doainance  cf  wavenumber  3,  7  and  11 

components  (Fig.  4)  . 

lhe  72-hcur  forecast  and  analysis  for  the  503rab  surface 
on  21  Feb  (Figs.  7a  and  b)  appear  guite  similar  under  a  cur¬ 
sory  inspection. 

The  trough  at  60-70W  has  been  under  forecast  by  the  model, 
but  overall  the  forecast  seems  to  be  reasonably  successful. 
Note  the  excellent  forecast  of  the  two  troughs  in  the 
Pacific  Ocean  (18QW  and  140H)  and  the  strong  ridge  over 
Furope. 

iher.  the  spectral  decomposition  is  applied  to  the  72 
hour  forecast  ar.d  its  verification  (Fig.  7c)  significant 
amplitude  differences  are  present.  specifically,  in  -.he 
planetary  wave  greup  (wavenumbers  1-3)  the  observed  atmos¬ 
phere  has  maintained,  although  to  a  lesser  degree,  the 
dominance  cf  wavenumber  3  whereas  the  model  has  shif-sd  the 
peak  amplitude  away  from  wavenumber  3  and  into  wavenumber  1. 
In  the  long  wave  group  (wavenumbers  5,6  7)  the  dominance  cf 
wavenumbers  6  and  7  has  shifted  to  the  sole  peak  in 
wavenumber  7  in  fceth  the  tcdel  and  the  analysis. 
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500MB  5FC  HEIGHT 
FEB  21 


NOGAES  accurately  forecast  the  wavenumber  amplitudes  in  this 
group.  The  medium  wave  group  (wavenumber  8-12  coapor.er.ts) 
initially  indicates  a  preference  for  wavenumber  11,  but 
after  72  hours  the  significance  has  shifted  to  the  lower 
wavenumbers  8  and  9.  This  shift  is  slightly  more  pronounced 
in  the  model  resulting  in  erroneously  high  forecast  ampli¬ 
tudes  for  wavenumbers  3  ard  9. 

The  Hovmoller  diagrams  for  the  planetary  wave  components 
of  the  forecast  and  the  analysis  (Figs.  3a  and  b)  clearly 
illustrate  these  differences. 
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Figure  8. 
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The  model  has  accurately  positioned  ridgss  at  120W  and  0 


but  it  has  underestimated  the  amplitude  or  the  former  by 


C 


Ths  rwe 


meters  and  overestimated  the  latter  by  50  meters, 
planetary  wave  treaghs  originating  at  18QH  and  40W  are  v~ak- 
ened  tco  rapidly  in  the  forecast  and  become  distorted  by  the 
spurious  formation  of  a  deep  trough  at  9QE.  The  trough 
originating  at  30E  also  weakened  toe  rapidly  to  the  point 
where  it  is  higher  than  the  mean  height  and  is  only  a  trough 
in  the  sense  that  it  is  between  two  ridgas.  The  ether  prom¬ 
inent  feature  is  the  forecast  trough  at  90E  which  has  no 
counterpart  in  the  observations.  Tha  seemingly  random 
changes  in  trough  and  ridge  intensity  form  a  pattern  when 
the  differences  in  component  amplitudes  ara  considered. 

The  NCGAPS  72-hour  forecast  shews  the  wavenumber  1  sig¬ 
nature  emerging  with  general  ridging  from  1 50 w  to  40E  and 
ganeral  troughir.g  from  40E  to  15i)W  while  tha  analysis 
depicts  mere  cf  a  wavenumber  2  and  3  signature.  Figs.  9a 
and  fc  compare  the  forecast  and  observed  wavenumbers  1  and 
individually. 

NOGAES  forecasts  wavenumber  1  to  be  70S  more  intense  than 
the  observation  while  underforecast ing  wavenumber  3  by  55%. 
Despite  the  presence  of  wavenumbers  2  and  3 ,  the  planetary 
wave  characteristics  in  the  72  hour  NOGAPS  forecast  ars 
largely  determined  by  the  overwhelming  dominance  of 
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Fiqure  9.  Wavenumber  1  ar.d  3  Components  ror  Casa  1.  a) 
wavenumber  1  Cciponen t  f or  Analysis  (Solid) 
Forecast  (Cashed)  for  OOZ  21  Feo.  b)  Same  a 
except  fcr  Wavenumber  3  Component. 


wavenumber  1.  Erroneously  lew  heights,  troughs  too  str cr.g 
and  ridges  toe  weak,  exist  in  the  region  of  60S  to  120W  and 
erroneously  high  heights,  troughs  too  weak  ar.d  ridges  tec 


(0  01 


strong,  exist  in  the  region  of  120fi  to  6QE.  Synoptic  sxam- 
ples  of  the  effect  of  the  reduction  of  the  wavenumber  3 
component  are  evident  in  the  ridge  at  12QW  and  the  trough  it 
60H.  Froi  150H  to  1 10H  the  model  under  forecasts  the  ridge 
in  this  area  and  smooths  cut  its  amplitude.  From  11GW  to 
60 a  the  model  o ver f o reca sts  the  height  of  the  eastern  side 
of  the  ridge  and  smooths  out  the  trouga  at  60S.  These 
amplitude  errors  are  also  evident  in  other  regions  of  the  72 
hour  forecast  (Figs.  7a  ard  b). 

The  long  save  group  verifies  the  wavenumber  5,6  and  7 
components.  The  general  agreement  between  the  individual 
wavenumber  amplitudes  in  this  group  (Fig.  7c)  through  the 
72 -hour  period  consequently  produces  a  general  agreement  in 
the  anplitudes  of  their  respective  Hovmoller  diagrams  (Figs. 
10a  and  b)  . 

This  diagram  shews  seven  troughs  ar.d  ridges  at  72  hours  sig¬ 
nifying  the  dcninance  of  wavenumber  7.  in  this  Grouping  a 
region  of  maximum  amplitude  (contour  intervals  are  15m)  is 
evident  between  18QW  and  30W  in  both  the  model  ar.d  the 
observation.  There  is  strong  agreement  between  the  forecast 
and  observed  long  waves  in  the  Hovmoller  diagrams  showing 
NOGAES  performing  well  at  this  scal9. 
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Figure  10.  long  Wave  Hovmcller  Diagrams  for  18-21  Feb.  al 
Forecast  and  t)  Analysis  with  Contour  Intervals 
of  30  Haters. 


The  interaction  of  waver.um bars  8  through  12  in  the 
medium  wave  group  creates  a  complex  pattern  in  the  analysis 


and  forecast  trough-ridge  diagrass.  The  Hovmoller  diagrams 
for  the  forecast  and  observation  (Figs.  11a  and  b)  both  show 
eleven  waves  initially  and  nine  waves  after  72  hours. 

The  anplitudes  (contour  intervals  are  10  meters)  ar.  d  phase 
speeds  can  still  be  compared,  but  tha  errors  are  acre  diffi¬ 
cult  to  trace.  The  model  shows  an  erroneous  less  of  ampli¬ 
tude,  particularly  evident  at  48  hours.  However,  at  72 
hours  the  increased  activity  in  the  region  from  140E  to  120W 
is  properly  forecast.  The  phase  speed  is  approximately  10 
degrees  tcc  fast  in  the  region  from  OOE  to  120E.  Con¬ 
versely,  the  model  has  done  extremely  well  in  forecasting 
the  speed  cf  the  wave  features  between  14QE  ar.d  150W.  Over¬ 
all,  the  mcdel  exhibits  reasonable  success  for  both  the 
amplitudes  and  phase  speeds  cf  the  medium  waves  at  72  hours. 

It  is  important  to  note  the  amplitude  of  the  medium  wavenum¬ 
bers  is  often  smaller  than  the  planetary  wavenumbers  by  as 
much  as  an  order  of  magnitude.  Also  as  wavenumbers  increase 
to  the  point  where  they  cannot  be  accurately  represented  by 
the  model's  horizontal  resolution,  the  spectral  decomposi¬ 
tion  is  less  accurate  as  veil. 

The  results  of  Case  1  indicate  NOGAPS  erroneously 
shifted  the  planetary  wavenumber  dominance  from  wave  3  and 

'j 
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model  and  the  observation  were  difficult.  Both  the  ampli¬ 
tudes  and  phase  speeds  cf  the  long  waves  were  accurately 
predicted  by  the  model.  NOGAPS  shewed  some  less  of  ampli¬ 
tude  ir.  the  medium  waves  and  shifted  peak  amplitude  from 
wavenumber  11  to  8  and  9.  The  phase  speed  errors  of  the 
medium  waves  as  a  group  were  small. 

4.2  Case  2  -  Cctcber  06  tc  11,  1982: 

This  period  provides  a  good  example  of  the  formation  of 
an  intense  ridge  (blocking  high)  over  the  western  U.  s.  and 
was  discussed  in  the  NOA  A,  Western  Regional  Attachment  of 
October, 1S82.  For  this  case  NOGAPS  horizontal  resolution  is. 
2.4  degrees  of  latitude  by  3.  degrees  of  longitude  interpo¬ 
lated  onto  a  2.5  by  2.5  degree  spherical  grid.  The  spectral 
decomposition  is  applied  tc  the  height  of  the  500mb  surface 
along  45N  (averaged  between  the  heights  given  at  42.  5N  and 
47  ,5N) . 

The  500mt  heights,  depicted  from  30N  to  60N,  for  06  Cct 
in  Fig.  12a  shows  a  fairly  zonal  pattern  from  the  east  coast 
of  Asia  to  the  Eastern  Atlantic. 

Nets  the  ordinate/abscissa  ratio  is  more  than  one.  This 
distortion  produces  a  visual  magnification  cf  the  amplitudes 
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Cct.  a)  500a t  Heights  with  Contour  Intervals 
120  Meters,  b )  Wavsr.u  mber  Component  Amplitudes 


Meters. 


of  each  wave.  This  scale  was  chosen  for  clarity  and 
maintained  throughout  this  and  the  remaining  cases.  Prom 


neat  ridges  are  located  at  20W  and  30E  and  the  flew  spli 


at  1 0H.  Spectral  decomposition  (Fig.  12b)  shews  wavenumber 
2,  5  and  7  components  are  dominant  initially.  The  zonal 
character  cf  this  case  is  shown  by  the  relatively  even  dis¬ 
tribution  cf  amplitudes  among  wavenumbers  1  through  8. 

Seventy-two  hours  later  (Fig.  13a)  a  deep  trough  has 
developed  ever  the  central  U.S.  and  prominent  ridges  have 
rapidly  developed  over  the  western  and  eastern  U.S.  coasts. 
The  model  (Fig.  1 3 1)  does  not  adequately  forecast  the 
strength  cf  these  features.  A  comparison  of  the  spectral 
decomposition  (Fig.  13c)  cf  the  analysis  and  the  forecast  at 
72  hours  shows  a  marked  increase  in  the  dominance  cf  the 
wavenumber  2  and  7  components  of  the  atmosphere  while  the 
model  incorrectly  smoothed  the  amplitude  of  wavenumber  2  tut 
correctly  maintained  ‘he  relative  strength  of  wavenumber  7. 

At  the  120  hour  mark  (11  Oct)  of  the  forecast  period 
both  the  eastern  and  western  U.S.  ridges  have  continued  tc 
strengthen  (Fic.  14a)  .  The  model  shews  an  ability  to  prop¬ 
erly  forecast  the  phase  speed  of  the  wave  features  such  as 
the  troughs  at  170W  and  40W  and  the  ridge  at  160E  but  con¬ 
tinues  to  underestimate  the  strength  of  the  ridges  (Fig. 
14b)  over  the  U.S.  Spectral  decomposition  of  the  observa¬ 
tion  shows  the  largest  amplitudes  are  in  wavenumbers  2,  3,  4 
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a)  500MB  SFC  HEIGHT 

NOGRPS  72  HR  FOR  OCT  09 


b)  500MB  SFC  HEIGHT 

OCT  09 


Figure  13 


and  5  all  having  similar  values  (Fig.  14c).  Meanwhile  -.he 
model  has  shifted  the  dominance  to  wavenumber  4  while  keep¬ 
ing  the  amplitude  of  wavenumbers  5  and  greater  relatively 
small.  The  largest  amplitude  changes  in  -he  model  during 
the  forecast  period  occurred  in  wavenumbers  2,  3  and  4  while 
the  atmosphere  additionally  showed  large  amplitude  changes 
in  wavenumber  5. 

The  s ignif icunce  of  the  amplitude  errors  is  readily  dis- 
cernatle  in  the  Hovmcller  plots.  Initially  bo-.h  cf  the 
planetary  wave  trough-ridge  diagrams  (Figs.  1 5a  and  h) 
exhibit  the  characteristics  of  a  large  wavenumber  2 
component. 

After  72  hours  the  observation  has  greatly  increased  the 
strength  cf  the  wavenumber  2  signature  while  the  forecast 
decreased  it  (Fig.  16a) .  Tha  model  has  successfully  fore¬ 
cast  the  strength  of  the  trough  at  160W  (-135m),  but  the 
strength  cf  the  other  trough  and  both  ridges  are  severely 
underestimated.  The  combined  dominance  of  wavenumbers  2  ar.d 
3  in  the  NCGAFS  forecast  has  introduced  a  spurious  ridge  at 
OH  and  trough  at  50E. 

Ey  120  hours  the  atmosphere  exhibits  maximum  amplitude 
in  the  planetary  wavenumbers  2  and  3  (Fig.  14c).  This  shift 
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Figure  14.  saae  as  Fig. 7  excapt  for  11  Oc 
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Figure  15.  same  as  Fig.  6  except  for  06-11  Oct. 


has  net  affected  the  amplitudes  cf  the  planetary  wave  fea¬ 
tures  as  much  as  it  has  sharply  tightened  the  spacing 


between  the  treugnes  ard  ridges  producing  acre  of  a 


wavenumber  3  signature.  The  model  has  forecast  the  ampli¬ 
tude  of  wavenumber  1  accurately  while  underestimating 
wavenumber  3  by  30%  (Fig.  16b). 

Thus  the  model  has  cnce  again  failed  to  predict  the 
shift  and  erroneously  smoothed  the  highest  amplitudes.  The 
model's  planetary  wave  Hcvmoller  diagram  shows  the  weaker 
amplitudes  at  120  hours  and,  significantly,  the  resultant 
smoother  pattern.  Specifically  the  ridges  at  105E  and  11QW 
are  too  lew  and  the  trough  at  45W  is  too  high.  Synoptically 
these  errors  are  evident  in  Fig.  14b  as  the  failure  of  the 
model  to  adequately  forecast  the  strength  of  the  ridge 
building  cn  the  west  coast  of  the  U.S.  and  the  trough  at 
40 w .  The  general  placement  cf  comparable  features  is  qccd 
especially  in  the  Pacific  region  near  the  dateline. 

The  long  wave  component  Hovmoller  diagrams  (Figs.  17a 
and  b)  show  NOGAPS  successfully  forecasting  the  long  wave 
features.  Overall  the  forecast  has  accurately  predicted  the 
general  regions  of  high  long  wavs  amplitudes,  120W  to  90E, 
ar.d  lew  long  wavs  amplitudes,  at  903  to  120».  Observed  long 
wave  amplitudes  are  more  intense  than  the  forecast  after  the 
72  heur  mark.  Comparison  of  model  and  observed  phase  speeds 
shows  NOGAES  to  be  fast,  particularly  in  the  high  amplitude 
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waves  between  120W  and  30H.  The  shift  of  the  atmosphere  to 
waverunber  6  signature  within  the  group  by  72  hours  was 
properly  forecast  by  the  model  and  wavenumber  5,  6,  ar.d  7 

components  have  less  than  10%  amplitude  error  (Fig.  13c). 

Ey  120  hours  the  aticsphere  shows  maximum  long  wavs 
amplitude  in  the  wavenumber  5  component  (Fig.  14c)  and  ths 
Hovmcller  plot  (Fig.  17a).  In  contrast  ths  model  (Fig.  17b) 
exhibits  sxall  long  wave  amplitudes  and  underestimates 
wavenumber  5  by  mors  than  50%  (Fig.  16c)  and  wavenumber  6  by 
30%  while  overestimating  wavenumber  7  by  60%.  Ths  resul*  of 
these  errcrs  is  a  general  smoothing  cf  all  ridges  and 


troughs.  Alsc  the  region  of  high  amplitude,  120W  to  90S, 
has  exhibited  slew  phase  speed  error  while  the  region  of  low 
amplitude,  90Z  tc  120s,  stewed  rapid  phase  speed. 

The  medium  wave  components  (Figs.  18a  and  t)  begin  the 
forecast  period  with  the  signature  of  wavenumber  8.  By  72 
hours  no  clear  wavenumber  preference  is  present  ir.  either 
the  forecast  or  the  atmosphere.  The  same  is  true  after  120 
hours.  The  initial  significant  amplitudes  of  these  waves 
are  lest  in  the  atmosphere  after  24  hours.  A  few,  such  as 
the  trough  originating  at  140W  and  the  ridge  at  120W,  main¬ 


tain  their  amplitude  throughout  most  of  the  five  day  period. 
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sysi:€ai  0»€I  the  western  O.S.  The  model  does  capture  the  new 
activity  at  the  five  day  mark.  An  overall  comparison  sug¬ 
gests  a  phase  speed  error  of  5-10  degrees  too  fast  where 


direct  comparisons  can  be  made.  Syr.opt icaily  this  differ¬ 
ence  is  net  visible  cn  the  scale  of  Pigs.  14a  and  b  tut 
would  show  up  as  a  slightly  rapid  movement  of  features  with 
a  wavelength  cf  30  to  45  degrees  over  the  120  hour  forecast 
period. 

Ihe  results  of  Case  2  indicate  KOGAPS  erroneously 
smoothed  a  shift  in  the  planetary  wave  activity  at  72  hours 
and  again  at  120  hours.  Planetary  wave  phase  speed  errors 
were  small  fer  comparable  features  and  negligible  in  the 
Pacific  region.  Although  long  wave  forecast  amplitudes  were 
smoothed,  NCGAES  correctly  produced  two  distinct  regions  cf 
contrasting  amplitudes.  Phase  speeds  were  slightly  too  fast 
in  the  high  amplitude  region  but  satisfactory  in  the  low 
amplitude  region.  The  amplitudes  cf  the  medium  wave  fea¬ 
tures  were  well  forecast,  but  showed  a  fast  bias  in  wave 
speed  and  several  key  features  were  lest  in  the  forecast. 

4.3  Case  3  -  Cctcber  22  tc  27,  1982: 

Curing  this  forecast  period  the  500mb  surface  changed 
from  a  predominantly  long  wave  flow  pattern  to  a  more  zonal 
situation.  In  particular  the  region  from  150P  to  140H 
underwent  the  breakdown  cf  a  blocking  ridge.  Again  the 
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NOGAES  horizontal  resolution  is  2.4  degrees  la-.  by  3 
degrees  long.  interpolated  onto  a  2.5  by  2.5  degree  spheri¬ 


cal  grid.  The  spectral  decomposition  is  applied  to  the 
heights  alone  45N  averaged  between  42. 5H  and  47. 5N. 

The  500mb  heights  from  30N  to  60H  for  22  Oct  (Fig.  19a) 
includes  prominent  troughs  at  63E,  125E  and  140W  with  lesser 
enoughs  at  75W  and  15W. 

Noteahle  ridges  exist  at  90E,  170E,  110Wr  454  and  20E.  The 
wavenumber  5  dominance  is  clearly  present  in  Fig.  19b  alone 
with  strong  amplitude  components  of  wavenumbers  3  and  4  due 
to  th sir  particularly  strong  trough-ridge  systems. 

After  72  hours  ampiitede  differences  are  already  appar¬ 
ent  between  the  500ab  patterns  of  the  model  and  the  analysis 
(Figs.  20a  and  b) .  The  model  has  underforecast  the  strength 
of  the  ridges  at  90S  and  35E  and  the  troughs  at  70E,  140E, 
140W  and  ICE.  The  only  obvious  NOGAES  overforecast  is  in 
the  ridge  at  1Q0W. 

The  individual  wavenumber  components  (Fig.  20c)  a4-  72 
hours  show  the  dominance  remains  with  4  and  5.  Large  ampli¬ 
tude  increases  have  occurred  in  wavenumbers  1,  1,  9  and  10 
and  large  decreases  at  2  and  3.  The  model  has  accurately 
predicted  most  of  the  changes  but  consistently 
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3*  C  130’  C  ISO*  C 
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SPECTRAL  DECOMPOSITION  OF  500HB  SFC 
ANALYSIS  22  OCT  AT  45N 


underestimates  the  strength  of  wavenumbers  6  through  10.  In 
terms  of  this  study's  wave  groupings  the  amplitudes  of  the 
planetary  waves  are  forecast  very  well,  the  long  waves  are 
too  weak  ty  10  to  20*  and  the  medium  waves  are  too  weak  by 
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30  tc  40* 


500MB  SFC  HEIGHT 
NOGRPS  72  HR  TOR  OCT  25 


At  the  120  hour  mark  both  phase  speed  ar.d  amplitude 
errors  are  evident  in  the  500mb  patterns  (Figs.  21a  ar.d  t) . 
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amplitudes  of  the  planetary  waves  are  too  strong  by  20  to 
30S,  the  long  waves  are  well  forecast  and  the  medium  waves 
are  tco  weak  by  30  to  40*. 


AMPLITUDE (METERS) 


500MB  Sr C  HEIGHT 
NOGRPS  120  HR  FOR  OCT  27 


vma 

IVJII 


imi 


wmvm a 


u*e  90*  c  iao*c  iso*c  ioo*o  iso*n  iao*a  »’**  6 


500MB  5FC  HEIGHT 
OCT  27 


■W1 

iK£fi^ii6gaq 

miisfjk'M 


U*  C  80*  C  120*  C  ISO'C  100*  M  ISO*  H  120*  M  90*  H  60*  H  X 


SPECTRAL  DECOMPOSITION  OF  500MB  SFC 
ANALTS1SIS0L1O)  RNO  120HR  FORECAST (DASHED) 
27  OCT  AT  45N 


3I.  l 


1.0  3.0  S.0  7.0  0.0  11.0  13.0  1S.0  17.0  19.0  21.0  21.0  20.0  27.0 

-  HAVE  NUMBER 


Figure  21.  Same  as  Fig.  7  except  for  27  Oct 


The  Kcvmcllsr  diagrams  for  the  planetary  wavs  components 
of  the  forecast  and  tbe  aralysis  (Figs.  22a  and  b)  both  show 
a  shift  from  the  sole  dominance  of  wavenumber  3  to  a  shared 
dominance  with  wavenumber  1  by  7 2  hours. 

The  forecast  phase  speed  error  of  the  model  is  negligible  at 
this  point.  The  wavenumber  3  features  are  forecast  too 
strong  as  a  result  cf  the  20%  overforecast  by  the  model. 
Examples  include  the  troughs  at  1103  and  150W  and  the  ridges 
at  QOS  and  30H.  At  120  hours  the  model  shifts  its  maximum 
activity  tc  wavenumber  3  while  the  atmosphere  shews  a  dis¬ 
tinct  wavenumber  1  signature.  How  the  NOGAPS  wavenumber  3 
features,  the  troughs  at  120E,  SOW  and  35E  and  the  ridges  at 
85W  and  35E,  are  much  toe  strong.  Fig.  23  graphically  pres¬ 
ents  the  1001  cverforecast  cf  the  wavenumber  3  component  by 
the  medal. 

The  positions  of  the  ridges  at  85a  and  35E  are  excellent, 
but  the  trough  at  12QE  has  been  erroneously  retrograded  15 
degrees.  Another  interesting  aspect  of  this  case  is  the 
retrogression  cf  the  planetary  wave  which  originally  was  at 
150H.  MOGAFS  handles  this  phase  shift  very  well. 

Hithin  the  long  wave  group  the  signature  of  wavenumber  5 
persists  through  72  hours  in  both  the  model  and  analysis 
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Pigure  22.  Saae  as  Fig.  6  except  for  22-27  Oct 


0.0  >30.0  -9041 


(Pigs.  24a  and  t)  .  K  general  ccaparison  of  observed  and 
forecast  long  oaves  shows  good  agreeaent  on  this  scale, 
particularly  in  phase  speeds.  The  aodel  underforecasts  the 
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Figure  23.  Wavenumber  3  component  for  Casa  3.  Same  as  Fig. 
9  except  Wavenumber  3  for  27  Oct. 


strength  cf  wavenumbers  6  and  7  producing  slightly  weak 
amplitudes  in  the  region  from  30S  to  120S. 

However/  by  120  hours  both  observed  and  forecast  cat- 
terns  show  a  cuel  dominance  in  wavenumbers  5  ar.d  6  with  the 
model  amplitude  a  little  too  weak  in  wavenumber  5  and  a  lit¬ 
tle  tco  strong  in  wavenumber  6.  In  the  last  24  hours  of  the 
forecast  period  the  atmosphere  introduced  a  sixth  trough  and 
ridge  by  a  general  retrograding  of  certain  existing  features 
(180H  to  4CE)  by  5  tc  1C  degrees  of  longitude.  The  model 
introduced  a  sixth  trough  and  ridge  by  a  rapid  eastward 
movement  of  certain  features  (25H  to  140E)  .  This  incorrect 
me*hcd  of  handling  the  dominance  shift  resulted  in  phase 
speed  errors  in  the  range  of  5  to  15  degrees  too  fast. 
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Figure  25.  Same  as  Fig.  11  except  for  22-27  Oct 


dominance  ty  120  hours.  3cth  diagrams  show  their  largest 
amplitudes  in  the  region  froa  80  E  to  120Sf.  Overall  the  mod¬ 
el’s  amplitudes  are  too  weak  as  a  result  of  consistently 


■-  * 


undetfcracasting  the  strength  of  the  medium  waves.  The 
phase  speeds  are  too  fast  by  approximately  10  degrees  in  the 
regicn  fros  2CE  to  160E. 

The  results  of  case  3  indicate  NOGAPS  erroneously 
shifted  the  planetary  wavenumber  dominance  into  wavenumber 
3.  The  error  was  large  enough  to  significantly  effect  the 
amplitudes  of  all  the  wavenumber  3  features.  NOGAPS  plane¬ 
tary  wave  phase  speed  was  excellent  at  72  hours  and  still 
good  by  120  hours.  The  forecast  amplitudes  of  the  long 
waves  were  only  slightly  smoothed  and  long  wave  phase  speeds 
were  excellent  at  72  hours  but  too  fast  by  120  hours.  The 
medium  wave  amplitudes  were  smoothed  and  phase  speeds  were 
toe  fast. 
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V.  CCNCiU^ICNS  aggs|^gQ|STI0NS  I2H  FUT03E 

In  this  study  one  three-day  and  two  five-day  NO  GAPS  500afc 
forecasts  were  spectrally  verified.  The  cases  were  chosen 
to  reflect  three  different  planetary  and  long  wave  situ¬ 
ations.  In  Case  1  the  5C2nb  pattern  is  characterized  by  a 
vigorous  three  wave  pattern.  In  Case  2  a  formation  of  a 
blocking  high  and  a  shift  to  a  low  index  flow  dominate  the 
500mt  pattern  while  Case  3  describes  the  breakdown  of  a 
block  and  a  shift  to  a  high  index  flow. 

Two  analysis  techniques  were  used  to  study  the  observed 
and  forecast  spectral  statistics.  The  amplitudes  of  each 
wavenumber  in  the  analysis  and  the  forecast  are  compared  at 
key  periods  during  the  forecast  and  Hovmollsr  diagrams  are 
prepared  fcr  key  wavenumber  groupings  following  Baumhefrsr 
and  Ecwney  (1578).  Kaver.umber  components  were  grouped  into 
planetary  (waves  1-3)  long  (waves  5-7)  and  medium  waves 
(waves  8-12)  and  displayed  on  trough-ridge  diagrams. 

In  all  three  cases  the  most  serious  errors  were  concen¬ 
trated  in  the  planetary  waves.  The  model  consistently 
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failed  to  correctly  predict  large  amplitude  charges  in  the 
individual  planetary  wavenumber  components.  NOGAPS  both 
predicted  large  amplitude  changes  which  did  not  occur  in  the 
observations  and  failed  to  predict  large  amplitude  changes 
which  did  cccur  in  the  analysis.  In  Case  1  the  model  fore¬ 
cast  an  erroneously  large  amplitude  increase  in  wavenumber  1 
while  in  Case  2  NOGAPS  smoothed  the  wavenumber  2  amplitude 
which  had  increased  dramatically  in  the  analysis.  Based  cn 
these  three  cases  a  systematic  error  can  not  be  established. 
The  erroneous  amplitude  shifts  caused  differences  in  the 
wavenumber  signatures  and  made  estimates  of  phase  speed 
errors  difficult.  Where  direct  comparisons  could  be  made, 
planetary  wave  phase  speeds  were  accurately  forecast. 

The  long  wave  group  (waves  5-7)  was  the  most  accurately 
forecast  by  NCGAFS.  Overall  amplitude  and  phase  speed  char¬ 
acteristics  were  well  predicted  although  amplitudes  were 
slightly  sicctfced  by  the  icdel  and  phase  speeds  tended  to  be 
slightly  toe  fast  in  the  regions  of  increased  long  wave 
activity  such  as  in  Case  2.  In  ether  regions  of  both  slow 
and  fast  movement  the  model  showed  a  high  degree  of  phase 
speed  accuracy. 
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The  medium  waves  (6-12)  verification  exhibited  some  sys¬ 
tematic  errors.  NOGAPS  tended  to  smooth  the  amplitudes  cf 
the  medium  waves  to  the  point  where  key  features  were  often 
lost.  Also  phase  speeds  of  the  medium  waves  were  generally 
too  fast  in  the  model.  Specifically,  examples  of  fast  phase 
speeds  in  Case  1  are  in  the  region  from  OE  to  120E,  in  Case 
2  from  1808  to  308  and  in  Case  3  from  2QE  to  160E. 

Several  suggestions  for  future  research  follow  from  this 
study.  First,  the  number  of  forecasts  spectrally  verified 
should  be  increased.  New  that  the  spectral  verification 
software  has  teen  prepared,  the  verification  could  be  done 
routinely.  The  goal  of  the  increased  number  of  cases  would 
be  tc  document  systematic  spectral  errors  in  nogaps. 

A  second  area  that  should  be  investigated  is  *he  rela¬ 
tionship  between  the  strength  of  the  circumpolar  vortex  and 
the  forecast  phase  speed  errors.  Variations  in  the 
wavenumber  zero  component  can  be  found  by  comparing  gradi¬ 
ents  cf  the  mean  height  of  the  forecast  to  the  observation. 
The  strength  cf  the  circumpolar  vortex  may  be  a  factor  in 
the  phase  speeds  of  the  shorter  synoptic  waves. 


A  third  suggestion  involves  the  question  cf  the  cptimui 
north-south  interval  for  verification.  This  study  used 
eight  degrees  of  latitude  for  Case  1  and  five  degress  for 
Cases  2  and  3.  Further  studies  should  examine  the  effects 
of  narrower  cr  wider  latitude  bands  in  the  analysis.  Simi¬ 
larly  the  effects  of  latitude  variations  or.  wavenumbers 
could  be  determined  by  using  data  from  different  latitudes 
(6  ON  ,  3 ON  ,  etc)  . 

This  study  illustrated  the  utility  of  diagnostic  verifi- 
caticn  techniques  for  study  of  operational  numerical  pre¬ 
diction  models.  Future  improvement  in  NQGAP3  forecast 
skill,  particularly  in  the  medium  range,  is  dependent  in 
part  cr.  the  continuing  study  of  the  model  to  isolate  its 
strength  and  weaknesses. 
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